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Abstract 

B.M. Mirza p] presented a solution of coupled Einstein-Maxwell equations for a slowly 
rotating neutron star; however his derivations had some errors and implicit assumptions 
that rendered the solution invalid. We point out the errors and present a mathematically 
consistent solution. The resulting solution is also physically consistent as it remains finite in 
the no rotation limit, whereas Mirza's solution diverges for zero rotation. 

Slowly rotating neutron stars were first investigated in 1967 as a slow rotation approxima- 
tion (0 < a < l,0(a 2 )) to the Kerr metric [2] and have astrophysical relevance as most of the 
observed pulsars are actually slowly rotating relative to the speed of light (one fifth for a millisec- 
ond pulsar) [3]. Charged particle dynamics around such stars is investigated by constructing 
the corresponding Einstein-Maxwell equations assuming the slow rotation approximation [4]. 
Later Mirza [T] solved the same model using an ansatz. The chosen ansatz was dimensionally 
inconsistent and yielded an unphysical answer for the fields which gives a divergent expression 
for the radiation emitted in the no rotation limit. We prove that his chosen ansatz, even after 
dimensional modification, gives a divergent result in the no rotation limit. We then suggest an 
ansatz that avoids divergences in any limit and is physically more meaningful. 
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The spacetime exterior to a slowly rotating neutron star is given by the slow rotation approx- 
imation to the Kerr metric: 



ds 2 = _ e mr) di 2 _ 2LO ( r )r 2 sin 2 Odtdip + r 2 sin 2 9dip 2 + e~ 2 ^dr 2 + r 2 d9 2 , (1) 

where 

e 2*(r) = (1 _2M } (2) 
r 

and 

/ \ _ d(p _ g ttp ^ 
dt Qiptf 

is the angular speed of a freely falling frame brought into rotation by frame dragging. Here and 
in what follows the Greek indices run as t,r,9, and (p respectively and we use gravitational units 
in which G = 1 = c. The general relativistic form of the Maxwell equations is: 

+ Fpi,a + F-ya,/3 = 0, (4) 



(y^ F ^) p = ^^j a , (5) 

where g is the determinant of the metric tensor g a p given by the Einstein equations and J a is 
4-vector current density. Here F a p is the generalized electromagnetic field tensor for an ideal 
fluid given by a unique tensor ial expression: 

F a p = u a Ep - upEa + r] a p jS u r B s , (6) 

where i] a ^s is the totally skew tensor, E a and B a are the electric and magnetic fields and 
u a is the unit velocity 4-vector [5]. In general J a is the sum of two terms corresponding to a 
convection and a conduction current: 



J a = eU a + au p F Pa (7) 

where e is the proper charge density, a is the conductivity of the fluid. For a zero angular 
momentum observer (ZAMO) u r and uq vanish, and using u a u a = — 1, the components of the 
4-velocity vector are: 

u « = e -*W(l } 0,0,w(r)) J u a = e* (r) (-l, 0,0,0). (8) 

The electromagnetic field outside the neutron star is now determined by eqs.(4) and (5). To 
solve this system of equations let us assume Mirza's ansatz modified to maintain dimensional 
consistency, for the electric field E: 
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E r (r, 9) = hE e (r, 9) = k 2 E^(r, 9) = R E {r)Q E {9), (9) 

where k\ and k% are some dimensional constants. For the magnetic field B, the corresponding 
ansatz will be 

B r (r, 9) = k 3 B e (r, 9) = k 4 B v (r, 9) = R B (r)Q B {9), (10) 

where k% and k 4 are some dimensional constants. Further, following Mirza, we take a constant 
angular speed of rotation, ab- 
solving eqs.(4) and (5) using eqs.(9) and (10), we get 

E r = k x E 6 = k 2 E v = R E {r)@ E {9) = ^— , (11) 

B r = k 3 B e = hB ip = R B (r)e B (9) = p—, (12) 

where A\ and A 2 are arbitrary constants. It can easily be seen that in the no rotation limit 
i.e u —> 0, the electric and magnetic fields become infinite, which cannot be true, because it 
yields infinite energy radiated by a non-rotating neutron star. Further, the Poynting vector 
S = E x B gives the momentum flux and hence the energy radiated e. It yields e oc — ^ where 
the proportionality factor depends upon ki.-.k^ and {r,9). If it is non zero, the radiated energy 
diverges as lo — > which is impossible, and if e = 0, the rotating star would not radiate and 
hence would not be a model for a pulsar. Infact, in general, for a rotating neutron star % oc 
p]. Hence Mirza's ansatz does not work even after correcting for dimensions. 

We present another ansatz that avoids the above mentioned impossibilities. 

E = (E r ,Q,E tp ),B = (0,B e ,Q), (13) 

where E and B has r and 9 dependence only. To solve eqs.(4) and (5) we shall use the following 
separation ansatz for electric and magnetic fields 

E r (r, 9) = k 5 E^(r, 9) = R E (r)Q E {9), (14) 
B e (r,9) = R B (r)Q B (9). (15) 
Solving eqs.(4) and (5) using eqs.(13), (14) and (15) we get 

E r = k 5 E^ = R E (r)e E (9) = 4^- t , (16) 

B e = R B (r)Q B (9) = ^-, (17) 
ut sin 9 
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where A3 and A4 are arbitrary constants. Hence, E and B remain finite in the no-rotation limit 
and the radiation vanishes for a non-rotating object, as required. 

There will be several other ansatz that would yield solutions to the system of equations but 
care must be taken to avoid non-physical solutions such as Mirza obtained. As an extension to 
the problem, one can solve the above system using the source terms i.e. j a / and deduce some 
physically interesting realizable results. 
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